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■ scattering in quantum wells for small deviations of the distribution function 
from equilibrium, within the framework of the Boltzmann equation. For a 
given temperature T and density n, a symmetric matrix needs to be evalu- 
ated only once, and henceforth it can be used to describe electron-electron 
scattering in any Boltzmann equation linear-response calculation for that par- 
ticular T and n. Using this method, we calculate the distribution function 
and mobility for electrons in a quantum-well, including full finite-temperature 
dynamic screening effects. We find that at some parameters which we inves- 



We describe a method for numerically incorporating electron-electron 
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tigated, electron-electron scattering reduces mobility by approximately 40%. 
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I. INTRODUCTION 



The effect of electron-electron (e — e) interactions on transport in bulk-like systems of 
various dimensionalities is still an area of active research to this day, both experimentallyBi 
and theoreticallyi'i. On the face of it, it would seem that in the case of doped parabolic 
band semiconductors where umklapp processes are negligible, e — e scattering should not 
affect the linear transport properties of bulk-like systems (purely quantum effects such as 
weak-localization corrections excepted) since an e — e scattering event conserves the total 
current in the system. Nevertheless, it has been appreciated for a long time that e — e 
scattering can affect the mobility of a system semiclassically by scattering carriers into or 
out of parts of the Brillouin zone which are strongly affected by the other available scattering 
mechanisms.!'! 

In Ref. ^, it was shown that the expression for mobility in the presence of quasi-elastic 
scatters takes on different forms in the limits of zero and infinitely strong e — e scattering. 
Given a quasi-elastic energy-dependent transport (i.e., calculated with the 1 — cos 6* term) 
scattering time T{e) due to other scattering processes in the system such as acoustic phonons 
or impurities, the mobilities of the system for the cases of zero and infinite e — e scattering 
rates, respectively, are given by 

e(r) 



m 



^^oo = I • (1) 



Here, 



where n is the carrier density of the system, eih) = fi^ k"^ / {2m) , fQ{e{k)) is the Fermi-Dirac 
distribution function, d is the dimensionality of the system, and we are assuming an isotropic 
parabolic band system. Clearly, for the case when that temperature T is small on the scale 
of the energy scale over which r^e) varies significantly, (r) ^ [{t^^)]^^- Conversely, in the 
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case where T is large on the scale over which T{e) varies, there can be significant differences 
in the calculated mobility using the two different methods. For example, we show below 
that in particular cases in GaAs quantum wells, /io//^oo ~ 2. Thus, for accurate theoretical 
determination of the mobility at the semiclassical level, it is important that e — e scattering 
effects are included. Furthermore, it has been shown0 that experiments measuring the drag 
rate between electron gases between two coupled quantum wells is sensitive of the exact 
details of the linear response distribution function / in each layer. Since / is strongly 
affected by e — e scattering in the intermediate temperature regime T ^ 0.5 Tp (where 
Tp is the Fermi temperature), it is important to include the effects of e — e scattering in 
calculations of the drag rate. 

In this paper, we demonstrate an efficient way of including e — e scattering in the cal- 
culation of for linear transport for two-dimensional cylindrically symmetric systems, within 
the semiclassical Boltzmann equation formalism. A similar calculation has been presented 
in three-dimensions.! Within this formalism, the Boltzmann equation for linear response 
can be solved exactly (within numerical accuracy). We have included full effects of finite- 
temperature dynamical screening, which automatically includes phenomena such as Landau 
damping and collective mode enhancements to scattering. The description of the e — e 
scattering formalism for the Boltzmann equation is given in sections || and |T|, and section 



rVl contains the results and discussion. Throughout this paper, we assume that bands are 
isotropic and parabohc. 



II. ELECTRON ELECTRON SCATTERING PROBABILITY 

The e — e scattering occurs in the presence of other conduction electrons, and hence 
the bare interparticle Coulomb interaction U{q) is screened. Furthermore, at the interme- 
diate temperatures in which we are interested, the energy transfer between the electrons 
in a scattering event is often a substantial fraction of the kinetic energy of the electrons, 
and hence the scattering matrix elements for e — e interactions should be calculated using 
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the dynamically scieened Coulomb interaction V{q,u!) = U{q)/e{q,uj), where e(g,c<j) is the 
dielectric function. In this paper, we use e{q,uj) given by the random phase approximation 
(RPA), which we evaluate using a method described previously by us,0 and we use the Born 
approximation for the scattering probability. 

The scattering probability w(fc'^cri, ^2(12; fci<7i7 fc2C"2) for a pair of electrons initially in 
states fciO"!, fc2a"2 to be scattered to k[cxi, ^2^2 depends on whether or not the electrons have 
the same or opposite spins. For electrons with the same spin, say t. 



1 2tx 

w{k + q],k' -q^;k^,k' t) = - y 1^(9, ^k+q ' ^fc) " " ^ + ^' ^k' q ~ ^k' 



(3) 



The fraction 1/2 in Eq. is due to double counting, since w{k + q ] , k' — q k ] , k' ]) 
and w{k' — q ] ,k + q ]\k ] ,k' ]) describe exactly the same process. For opposite spins, 

w{k + qtk'-ql;ktk' i) = Sf^^^ - 6^)]'. (4) 

There is an equal probability that an electron scatters off another electron with equal or 
opposite spin, so one can sum over the Eqs. (Q) and (|) to obtain an "average" scattering 
probability! 

27r r 

w{k + q,k' - q;k,k') = ^1 1^(9' ^fc+g " ^fc) 1^ 

- lMV{q, Bk^q - Sk) V*ik' -q-k, Sj^,^ - e^)]]. (5) 

The first and second terms in Eq. referred to as the direct and exchange terms, re- 
specitively. 

In practice, the exchange term often makes calculations considerably more complicated 
and is usually ignored. The physical grounds for doing so are as follows. First, e — e collisions 
are usually dominated by small q scattering (because V{q, e) falls quickly with q for finite 
e). The direct [exchange] term has the form e:)^ \V {q, e)V {q' , e') with q' ^ g], which 
implies that the direct term usually dominates over the exchange. Also, the sign of exchange 
term can sometimes be negative, which leads to cancellation of this term within collision 



integral. The effect of the exchange term was studied in a 3D system with scatically screened 
interaction.! There it was found that the exchange term was significant for na^ = 1 but not 
for na^ = 0.1. In the calculations that follows we have used na"^ = 0.15, and since we 
furthermore include dynamical screening, which leads to a peaked interaction at small q, we 
can assume that the direct interaction dominates in our case. 

We write down the formal expressions for the electron-electron scattering operator both 
including and excluding the exchange term. However, in the actual numerical evaluation, 
we ignore the exchange interaction. 

III. ELECTRON ELECTRON SCATTERING OPERATOR 

The Boltzmann equation for electrons in uniform electric field E producing a force F = 
{-e)E is 

\ / e—e \ / p,i 

where the subscripts e — e and p,i are for scattering due to electron-electron interactions 
and the phonon+impurity interactions, respectively. 

We define the function \E'(fc), related to the deviation of the distribution function from 
equilibrium, as 

f(k)-fo{k)^fo{k)[l-fo{k)]^{k). (7) 
This function can be written in terms of a sum of angular components 

^{k,e)=J2^/j4k) cos(n^). (8) 

n 

where 6 is the angle from an axis of symmetry (here, the direction of the electric field). By 
the assumption of cylindrical symmetry of the system, the scattering terms in the Boltzmann 
equation do not mix different cos{n6) components, and one can isolate and concentrate on 
the cos^ component, ipi{k), since this is the one that affects the current and hence the 
mobility. 
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The cos^ component of the hnearized electron-electron collision operator (i.e., neglecting 
higher powers in '^{k)), which we denote Ie-e['ipi\, is0 

Ie-e[i/JiXf.^p]{k) = -2 y J -^y^ w{k + q,k' - q, k, k') x 

f{k)f{k')\l-f\k + q)\ \l-f{k'-q 



^ \^k + ^k' - ^k+q - q 
{^i(A;)xy^ _p + ipi{k')xf^' p 

- M\k + Q\)xk+k,F - V^id^' - (9) 

Here, ^ is the cosine of the angle between k and fc'. The goal is to write the operator / 
in the form 

I[ipiXj^ p]ik)=X]^ p dp p K{k,p) ipi{p). (10) 

The kernal K{k,p) is symmetric, from detailed balance,0 and the extra factor of p in the in- 
tegral comes from phase-space. Thus, in order to incorporate electron-electron scattering for 
a particular density and temperature into a calculation, one need only generate K{k,p) once 
and store it; it can then be used for all calculations involving electron-electron scattering at 
that density and temperature. 

The four tp^s in Eq. (P) give four terms, each of which give a contribution to the kernel, 
K = Ki+ K2 + + i^4. In the following subsections, we explicitly write down the form of 
each of these kernels. 

A. First term, involving ipi{k) 

The '4'i{k)xj^ p can be factored out, and we obtain 

K,{k,p) = -5{k-p)k-'Uk) J [1-Mk + q)] 

r dk' 

2 j J^foik') [1 - /o(fc' - q)Mk + q, k'; k, k' + q) 
^ (^k'+q - ^k' - i^k+q - ^k}) ■ (11) 



In the event where the exchange interaction can be neglected, one obtains, as in Ref. |^ 
(we denote the scattering integral which neglects the exchange interaction with an asterisk) 

rf{k) f dq „n., 27r, 



KKk,p) = -S{k -vi^jj^A^- + q)] -^\V{q, Sk^q - Sk) 



Tik J (27r)^ 



(12) 



where ub^s) = [exp{Pe) — 1] ^ is the Bose function, and x{Q:^) is the RPA polar izability. 



B. Second term, involving ipi{k') 



Since cos(0 + 9') = cos{6) cos{6') — sin(^) sin(^') and the sin terms vanish from symmetry 
considerations, we can write xj^' p = F^k' k' Then, the second kernel is 

K2{k,p) = —2 — 7—77: — / dOu ,^ cosO 



(13) 



(27r)2 Jo ""^'P ^^^^k,p 
dq 



^w{k + q,p~ q- fc,p)[l -f{k + q)][l -f{p- q)] 



(2vr) 

^i^k + ^p-^k+q-^k-q)- 
The q integration can be evaluated by the change of variables 

q = Q + Ak/2, 
Ak = p — k, 
k = {k + p)/2. 



(14) 



Then, the 5-function in Eq. (|13[) becomes 



6{ef^ + ep- Ef^ q - 6p q) = 6 



m 



Ak' 



which gives 



K2{k,p) = ^-^^ dd^p cos%p 



2-K 



w{k + Q,k - Q;k + Ak, k - Ak) 



where (j) is the angle between Q and k. 



(15) 



(16) 



C. Third term, involving ipi^k + q\) 



Using Xf^^q F ~ F^k k+q ™ integrand and letting p — k + q, gives 

i^3(fe,p)- '^°^^|y^^V o'^^^P.fe / ^J'{k'+p-k)[l-f{k')] 

w{p, k'; k,k' + p- k)6{sp + Sf^' -e^- ^k'+p-k^' (^^) 

The 5-function 5{h? {p—k)-{k' —k) / m) reduces the dk' dimensional integral to one dimension. 
If one neglects exchange, then as with the first term the dk' integral can be done, giving 

1 [f\ep) - fiej,)] 



2n^Asmh'[{ep-e]^)/{2kBT)] 

rir 271 

I ^^k,p^^^^k,p X'^^*^ -^^^k- ^P)\^ Im[x(p -k,Sf^- Sp)] (18) 



D. Fourth term, involving ijji{\k' — q\) 

The kernel is 

2/0 (A;) /-stt 
^4(A:,p) = d0pkcos{epk) 

y-^w{k + q,k'- q; k, k')[l - f{k + q)] f{p + q) 
6 {ej^ + 6p+q - Sf^^q - sp) } (19) 

The term in the 5-function goes as 

^p+q-£p + £k-^k+q = —iliP-^)}^ (20) 

which reduces the q integration down to one dimension. 

In fact, the kernels for higher order components are very similar to the ones given above. 
For an angular variation proportional to cos n9, the Ki term is identical for all n, whereas 
with K2, Ks and K4 one simply replaces cos^ with cosn^ in the 9 integration. 

We have shown that one can calculate the matrix K{k,p) which gives the electron- 
electron scattering term for small deviations from the equilibrium. Once the matrix K[k,p) 
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has been calculated, one simply needs to iterate the equation for ipi^k) until convergence is 
obtained. 

In order to calculate ipiik) for the case when the e — e scattering rate dominates, it 
is often useful to use the fact that electron-electron scattering leaves a drifted Fermi-Dirac 
distribution invariant.0 Thus, for the case of elastic or quasi-elastic collisions, one can define 

Mk)=Mk)-^7ik) (21) 

where iljf^{k) = hkvd/ {ksT) is for a drifted Fermi-Dirac distribution. Any Vd can be used. 
In our case, because we cut off the matrix K{k,p) at a fcmax, which implicitly sets ipi^k > 
kmax) = 0, we chose Vd so that ^/'i(fcmax) = so that the distribution function is continuous 

at fcmax- 

We write the linearized e — e scattering term as 

I^_^[ij,{k)] = + J[Mk)] (22) 

where the first term on the right hand side corresponds to the diagonal Ki term and J 
corresponds to K2 + + K^. The Boltzmann equation for ipi{k) in the case when the other 
scattering mechanisms are quasielastic (which might include acoustic phonon scattering, 
which generally involves very small energy electron loss) is 

„i,jah\ h(k)\i-h(k)]{Mk) + <i,Y^(k)) Uk) ^„j,i,, 

eEv(k)\^-^= ^ + <2^' 

where T~i^{k) is the quasi-elastic scattering rate. This implies that one must iterate the 
equation 

~ . _ ^/o(fc)[i - /o(fc)] + n^i] - um - u{k)]^^^{k)r:,\k) 

"^'^ ^ ~ urn - fo{k)]r-,\k) + r-_\ik) ^'^^ 

to find ipi{k). 



IV. RESULTS AND DISCUSSION 



We study the case of electrons confined in a 100 A wide square GaAs quantum well with 
infinite barriers. We assume that the there is a 5-doping layer of (uncorrelated) charged 
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impurities, equal in density to that of the electrons in the well, situated a distance d away 
from the center of the well. We included three scattering mechanisms: e— e, charged-impurity 
and acousitic-phonon scattering, and we approximated the acousitic-phonon scattering as 
being elastic. 

We calculated the matrix in the form on 200 by 200 grid-points from k = to k = 5k p, 
and we used spline routines to interpolate between the grid points. The K3{k, p) and K^^ik, p) 
diverge logarithmically as — ^ p, which complicates the splining procedure, but we got 
around this problem by splining K3^4{k,p)/ ln{\k —p\), which is a smooth function. 

In Fig. 1, we show the deviation function ipiik) / \fk for a fixed density n = 1.5 x 10^^ cm~^ 
and temperature T = 30 K, but with several different distances d of the ionized impurities 
from the center of the quantum well. Note that when the distribution is a drifted Fermi- 
Dirac function, ipi(k)/\^ = constant. Thus, as the impurities are moved further away, the 
impurity scattering becomes weaker and the e — e scattering starts to dominateB and drive 
the distribution function closer to a drifted Fermi-Dirac function. The inset shows ipi{k) 
calculated both including and excluding electron-electron scattering for d = 100 A, which 
shows more clearly the effect of e — e on ipi{k). While transport experiments in a single 
layer are not particularly dependent on the details of the shape of ipi{k), it has been shownB 
that drag experiments in coupled quantum wells are quite sensitive to the details of ipii^k). 
In particular, when V'i(^) rises faster than (which implies that there are more carriers in 
the high energy region than for a drifted Fermi-Dirac distribution), the drag rate increases 
because high energy particles give a larger contribution to the overall drag rate, and there 
is greater opportunity for coupling to the plasmons of the system, which also enhances the 
drag rate. Therefore, for the purpose of calculating the drag rate in coupled quantum wells 
in intermediate temperatures, it is crucial to calculate the actual form of ipiik) accurately, 
including all salient scattering mechanisms. 

Fig. 2 shows the mobility /i as a function of ionized impurity distance d from the center 
of the quantum well. Also shown are the mobilities /xq = e{T) /m and /ioo = e/((r~^)m), for 
the limits of zero and infinite electron-electron scattering, respectively. The /io is generally 
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larger than /ioo because e — e scattering tends to scatter "runaway" electrons with large 
velocities (where the impurity scattering rate is small) back into lower velocity states. The 
inset shows that fio for this case can be almost twice fioo- 

As d becomes larger, the e — e scattering dominates over all other scattering mechanisms 
and fi fioo- Conversely, for small d, the impurity scattering is relatively large compared 
to the electron scattering, and /i is closer to the /iq than /ioo- The crossover from /iq to 
/ioo is shown with the open squares in the inset of Fig. 2. The crossover occurs when the 
impurity scattering and electron-electron mean free paths become equivalent. The transport 
scattering rate, for impurity scattering is given by r^^^ = ^/ {l^fn) = 2.6 x lO"*^^ cm^V^^s^^/ /i. 
The electron-electron scattering in two-dimensional systems is approximately given 



' e—e ^ 



n 2 



In 



E, 



+ In 



,Qtf 



Pf 



+ 1 



(25) 



For this system, this is on the order of 10^^ s^^. Thus, the crossover point, which should 
occur when these two are equal, is given by /i ~ 3 x lO^cm^V^^s^^. An inspection of Fig. 
2 also shows this to hold. Chabasseur-Molyneux et a/.i have experimentally found this 
crossover in GaAs/AlGaAs heterojunctions. Finally, at d = 150 A, fi ^ 0.6/io, implying the 
e — e scattering has caused a substantial reduction in the mobility. 

To summarize, in this paper we have described a method of including electron-electron 
scattering, including full finite-temperature dynamical screening, exactly in the Boltzmann 
equation, for small deviations of the distribution function from equilibrium. Using this 
method to calculate the distribution function and mobilities for electrons in a GaAs quantum 
well, we find a well-defined crossover from fio to /ioo (which can be significantly different from 
each other) when the e—e and impurity scattering mean free paths are equivalent. For certain 
parameters studied, e — e is responsible for reduction in the mobihty of up to 40%. 
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FIGURES 

FIG. 1. Deviation function il)i{k)k'~^/'^ (normalized to 1 at A; = 0) for electrons responding 
to a weak static electric field in a GaAs quantum well, width lOOA, density n = 1.5 x 10^^ cm~^ 
and temperature T = 30 K, for distances d = 50,150,250 and 350 A of the charged impurity 
layer from the center of the well. The collision term includes screened impurity, acoustic phonon 
and electron-electron scattering terms. The further the impurities are moved away from well, the 
more dominant the electron-electron scattering becomes and the distribution tends to a drifted 
Fermi-Dirac (a pure drifted Fermi-Dirac is a straight horizontal line). The solid (dashed) curve 
in the inset shows the deviation function (in arbitrary units) for d = 100 A calculated including 
(excluding) electron-electron scattering. 

FIG. 2. Mobility as a function of impurity distance from the center of the well. Parameters as 
in Fig. ||. Crosses are for ^ calculated from the Boltzmann equation, solid line is for /xq (no e — e) 
and dashed line is for /Xqo (infinite e — e). Inset: (i) (/i — ^oo)/(/^o — /^oo) as a function of impurity 
distance from center of well (solid line). As the mobility passes through 3 x 10'^ cm^y~^s~^, 
where e — e and impurity scattering rates are approximately comparable, the mobility crosses over 
from to /Uoo- (ii) The ratio fJ-oo/fJ-o (dashed line). This ratio can be as small as ~ 0.5. 
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